Vortex Lattice Melting and H c2 in underdoped YBa 2 Cu 3 
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Vortices in a type-II superconductor form a lattice structure that melts when the thermal displace- 
ment of the vortices is an appreciable fraction of the distance between vortices. In an anisotropic 
high — T c superconductor, such as YBa2Cu30 H , the magnetic field value where this melting occurs 
can be much lower than the mean-field critical field H C 2- We examine this melting transition in 
YBa2Cu30 H with oxygen content y from 6.45 to 6.92, and fit the data to a theory of vortex-lattice 
melting. The quality of the fits indicates that the transition to a resistive state is indeed the vortex 
lattice melting transition, with the shape of the melting curves being consistent with the known 
change in penetration depth anisotropy from underdoped to optimally doped YBa2Cu30j,. From 
the fits we extract H C 2(T = 0) as a function of hole doping. The unusual doping dependence of 
H C 2(T = 0) points to some form of electronic order competing with superconductivity around 0.12 
hole doping. 



I. INTRODUCTION 

Cuprate high-T c superconductors are of great interest 
not only because of their high transition temperatures, 
but also because strong-correlation physics gives rise to 
peculiar normal-state properties. Ironically, however, the 
strength of the superconductivity in these high— T c mate- 
rials is what interferes with measurement of the normal 
state properties at low temperature. Applying high mag- 
netic fields can overcome this and has led to the discovery 
of a small Fermi surface in underdoped YBa2Cu3 0, / via 
quantum oscillation measurements in pulsed fields. 1 This 
discovery prompted a large experimental survey of the 
transport and thermodynamic properties of YBa2Cu30 y 
in high fields. The questions remain as to whether the 
high fields are revealing the normal-state properties of 
YBa2Cu30y, or are instead exposing a qualitatively dif- 
ferent field-induced ground state, or whether one might 
still be in a regime dominated by superconducting pairing 
and the presence of vortices. 

The idea of high magnetic fields revealing the normal 
state in cuprate superconductors is a contentious one, 
in part because the phase diagram of the cuprates dif- 
fers qualitatively from that of conventional type-II su- 
perconductors. Owing to the short coherence length, low 
superfluid phase stiffness, and strong anisotropy, fluctu- 
ations play a dominant role in the phase diagram. There 
is evidence for 3D-XY critical fluctuations below and 
above T c . 2 4 Previous transport measurements on several 
cuprate compounds have shown that reaching the resis- 
tive state requires very high magnetic fields, and that the 
onset of resistivity as a function of field and temperature 
does not follow the conventional H C 2 curve derived from 
Ginzburg-Landau theory, as it does in more conventional 
type-II superconductors. 5 6 Instead, as is expected for a 
superconductor governed by strong thermal fluctuations, 



a vortex melting transition occurs, 5 7 8 with an extensive 
crossover regime to the normal state. Some Nernst effect 
experiments have been taken as evidence for the presence 
of superconducting pairing far above T c , even in strong 
magnetic fields. 9 With this in mind, it is important to 
consider at which field scale is superconductivity com- 
pletely suppressed and the normal state recovered, es- 
pecially with regard to quantum oscillation experiments 
which are purported to probe the "normal state" Fermi 
surface. 



II. THEORY 

The thermodynamic critical field H c is the field at 
which superconductivity is destroyed in a type-I super- 
conductor, and is directly related to the condensation 
energy of the superconducting ground state. In a type- 
II superconductor the magnetic field can penetrate the 
sample at a field lower than H r . At this field, H cl , the 
magnetic field penetrates the superconductor in the form 
of vortices, with each vortex being supercurrent running 
around a normal state core and containing a quantum of 
magnetic flux. The cores of these vortices, whose size is 
of order the superconducting coherence length £o, are in 
the normal state; outside of the vortex cores, the strength 
of the magnetic field decays over the length scale of the 
penetration depth A which, for strongly type-II super- 
conductors such as the cuprates, is much larger than 
the coherence length. These vortices can form a two- 
dimensional lattice perpendicular to the applied field (a 
"vortex lattice"), 10 and the lattice spacing shrinks in size 
as the magnetic field is increased. As long as the vor- 
tices remain pinned, the zero-resistance property is main- 
tained in the material. When the vortex cores overlap at 
a second field scale H C 2, superconductivity is destroyed. 
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In an isotropic, low-T c type-II superconductor, such as 
Nb 3 Sn, resistivity sets in at H C 2 and the diamagnetic sig- 
nal of superconductivity completely disappears. In terms 
of the mean-field Ginzburg-Landau coherence length 
this field scale is 



HoH c2 (T = 0) = 



_$o_ 



(1) 



where $o is the flux-quantum in SI units (H C 2 (T — 0) will 
henceforth be £f c2 (0)]P. 

The situation is more complicated in high-T c materials, 
where the vortex lattice can melt into a vortex liquid well 
below H C 2 . The Lindemann criterion for melting requires 
the thermal displacement of a lattice to be some frac- 
tion (defined cl) of the average lattice constant. Using 
the Lindemann criterion for a vortex lattice, Houghton 
et al. 12 have shown that, because of the large anisotropy 
in the cuprates, the vortex lattice in a strongly type-II 
superconductor with a high T c can melt at field values 
B m well below H C 2 for intermediate temperatures (away 
from K and T c ). 13 In these materials, H C 2 represents 
a crossover from a vortex-liquid to the normal state. In 
the traditional picture the melting field B m and /j,qH C 2 
are equal at zero temperature, since there are no thermal 
fluctuations at zero temperature to melt the vortex lat- 
tice. The presence of strong quantum fluctuations could 
result in a vortex liquid persisting down to zero temper- 
ature. However, in order to compare our experimental 
data with the theory of vortex lattice melting, we use the 
assumption made by Houghton et al. 12 , Blatter et al. 13 , 
and others that B m (0) = /io#c2(0)- 

Using the notation of Blatter et al. 13 , the melting tran- 
sition field B m is given implicitly by 



1 - M*) \f^ zr t 



4(>/2-l) 



2ttc| 
VGi' 



(2) 



The reduced field variable is b rn = B m / /xo-ff C 2j and t — 
T/T c is the reduced temperature. 

The Ginzburg number Gi, on the right hand side of 



Equation 2 is given by 
1 / k B T c j 



Gi = 



2 ^(rt(T=0)) 2 « 



(3) 



(9.225 x 10 8 [Wb^K- 1 ] x ^H c2 (Q)T c \ ab \ c y 



where 7 is the anisotropy ratio 7 



(4) 



and the defini- 



tion H c2 (0) 



A c are the penetration cfepths parallel and perpendicular 
to the a-6-plane at zero temperature). As emphasized by 
Blatter et al. 13 , this Ginzburg number should be thought 
of as a useful collection of parameters, and not as a num- 
ber describing the width of fluctuations around T c as it 
is in more three-dimensional superconductors. The Lin- 
demann number cl appearing on the right hand side of 



has been used (A & and 



Equation 2 represents the fraction of the vortex lattice 

that the thermal displacement 



parameter, a v = 

must reach in order for the vortex lattice to melt. 7 12 13 
Attempts have been made to calculate cl, with values be- 
tween 0.2 and 0.4 obtained for the cuprates, depending 
on the specific model (see Blatter et al. 13 for a review), 
but cl is probably better left as a fit parameter. 14 



III. EXPERIMENT 



Figure 1 shows a typical set of c-axis resistivity curves 
up to 60 tesla, from 1.5 to 200 K for YBa 2 Cu 3 6 . 5 9- We 
define the resistive vortex-melting transition as the mag- 
netic field where the resistance is l/100 t/l of its value at 
60 tesla. The definition of B m from resistivity curves is 
somewhat uncertain because of the width of the resistive 
transition (see upper panel of Figure 2). An alterna- 
tive definition would be the intersection of a line tangent 
to the steepest part of the resistive transition with the 
temperature axis. This would lead to small offsets (one 
tesla at most) in B m , but would not otherwise affect the 
conclusions of this paper. However, it is important that 
a consistent definition across different doping levels be 
used. 

The upper panel of |Figure 2| shows the vortex lattice 
melting transition from 1.2 K up to T c for YBa 2 Cu 3 0g.5g, 
one of the underdoped samples in which the melting tran- 
sition is accessible even at low temperatures. The con- 
cave upwards shape is characteristic of a vortex melting 
transition, as seen before in YBa 2 Cu 3 O y and in other 
cuprates, 6 14 and differs qualitatively from the concave 
downwards curvature of i? c2 (T) in conventional super- 
conductors. This form has been observed in a num- 
ber of cuprates, 6 15 16 but a systematic comparison to 



Equation 2 across the underdoped regime of the cuprates 
has not been performed. Here we present data for 
YBa 2 Cu 3 O y from oxygen content 6.45 to 6.92, with T c s 
ranging from 44.5 to 93.5 K, and identify trends that arise 
as a function of doping. Characteristic curves for several 



other dopings are shown in Figure 3 all with an upwards 
curvature, although that shape becomes less pronounced 
for the higher T c samples. 

|Equation 2| can be expanded about T c and solved for 
B m as shown in Blatter et al. 13 , but if the full tempera- 
ture range from 1.5 K to T c is to be used then it is more 
accurate to fit to the full implicit expression for B m . The 
use of |Equation 2| requires both the in and out-of-plane 
zero-temperature penetration depths, as well as the T c : 
these values are also listed along with th e hole doping 
(estimated using Liang et al. 17 ) in Table I The in-plane 
penetration depth values, A a t, come from electron-spin 
resonance (ESR) measurements 18 and from muon-spin 
rotation experiments, 19 both performed on comparable 
YBa 2 Cu 3 0j, crystals grown at UBC. In the case of the 
ESR values, the geometric mean of A and At, was taken. 
Out-of-plane penetration depth values, A c , come from 
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FIG. 1. The c-axis resistance of YBa2Cu306.59 as a function 
of magnetic field, from 1.5 to 200 K. The onset of resistivity as 
field is increased marks the vortex lattice melting transition. 
At low temperatures, quantum oscillations are seen above this 
melting field. 
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FIG. 2. Top panel: A magnified plot of the resistive tran- 
sitions shown above in |Figure l) The red dots are where 
the resistance is 1/100 4 of its value at 60 Telsa (extrapo- 
lated for high temperatures where resistance was not mea- 
sured to the highest fields). Bottom panel: The same data 
points as highlighted in red in the top panel, now plotted as 
a function of temperature. The black line is a fit is to |Equa^| 
|tion 2| using the known parameters given in |Table"l} and gives 
VoH c2 (0) = 28 ± 0.3 T, and c L = 0.37. 
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FIG. 3. The vortex lattice melting transition for several dif- 
ferent oxygen concentrations. The temperature axis has been 
scaled by T c , and the lines are best-fit lines to |Equation 2| 
All of the data points were aquired in the same manner as 
described in the caption of |Figure 2| 



infrared reflectance measurements, 20 also performed on 
UBC crystals. Interpolated values for the penetration 
depth were used when the exact doping values were not 
available. The penetration depth values and the interpo- 



lation are shown in the upper panel of Figure 4 

With A c , A a fc, and T c experimentally determined, the 
data at each doping can be fit using only two parameters: 
cl and H c2 (0). The fits in the lower panel of Figure 2 
and in |Figure~3| clearly show that three-dimensional vor- 
tex melting describes the in-field resistive transition in 
YBa2Cu30 y from y = 6.45 to 6.92. The penetration 
depth anisotropy, 7 = changes from ~ 50 at 6.45 
to ~ 16 at 6.92: this results in decreased curvature of 
the melting line as oxygen content (and hole doping) in- 
creases. This is the same behaviour seen in several dif- 
ferent cuprates of varying anisotropy, reported in Ando 
et al. e . The cl and H C 2(0) values extracted this way 
are given in |TableT] for all of the dopings measured. The 
fact that the Lindemann number remains relatively con- 
stant as a function of doping means that the shape of the 
melting curve is determined primarily by the penetration 
depths, which are becoming less anisotropic as hole dop- 
ing increases. The Lindemann number and the penetra- 
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Equation 2 



tion depths appear only as the ratio t — y 
and we plot this ratio in the lower panel of|Figure~4| The 

c 2 . . . 

increase of A ^ with hole doping is what is controlling 
the changing curvature as a function of doping. With 
this parameter setting the shape, H C 2(0) corresponds to 
the T = intercept of the melting curve. 



It should be emphasized that in Equation 2 



lim B m (T) 

T-+0 



(5) 



and so the values of H c2 (0) derived from fits to Equa- 
tion 2| are determined mostly by the zero-temperature 
intercept of the data for B m vs. T, and are essentially 
independent of the penetration depth values chosen. The 
penet ration depth s and the Lindemann number always 



enter 



Equation 2 



as the ratio 



and so errors in the 



penetration depth values (which arise because we use in- 
terpolated values from the upper panel of Figure 4 ) are 
absorbed into the fit value of cr, . 



In using jEq uation 2[ we have ignored the possibility 
that the onset of finite resistivity is due to the lattice 
dcpinning, and not actually melting. This assumption 
is probably justified, as the dcpinning transition is dis- 
tinct from the melting transition in YE^CusOj, only for 
temperatures very close to T c and in sa mples with ex- 
tremely low disorder. The fits shown in Figure 2 and 
|Figure 3l show data up to near T c when available, but only 
data at temperatures less than 0.8 x T c were used in the 
fits (which also avoids any possible effects of XY-critical 
phenomena near T C )J Additionally, all of the samples in 
this study (except for possibly the YBa 2 Cu3 06.92 sam- 
ple) have more disorder than the YBa2Cu30g.95 sample 
used in Liang et al. 5 . This is because the ortho-II, III, 
V, and VIII states are not perfectly ordered, 21 and have 
more disorder than ortho-I ordered YBa2Cu306.g5, which 
is close to stoichiometry. This disorder pushes the depin- 
ning transition closer to T c . 

The extracted values for -ff C 2(0) are plotted with the 



phase diagram of YBa2Cu 3 !/ in Figure 5 and show an 
anomaly around 0.12 hole doping. The solid blue line in 
Figure 5] is the function 



i - TcAT 



82.6(p-0.16) 5 



(6) 



where T™ ax is the maximum T c of the material (equal to 
94.3 K for YBa2Cu30j,). 17 This function has been found 
to describe T c as a function of p in the cuprates, except 
for the suppression of T c around l/8 th hole doping. 17 
The green circles in |Figure 5|are the abso lute difference 
between the actual T c and |Equation~6{ and the sup- 
pression of T c is clearly correlated with a suppression 
of H C 2 (0) . Suppression of the melting transition in this 
region was reported for a few different doping levels in 
LeBoeuf et al. ^. 
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FIG. 4. Upper panel: The in and out-of-plane penetration 
depths of YBa2Cu30 !/ , as measured by muon-spin rotation 
(X a t at 0.11 holes), electron-spin resonance (the other X a t 
points), and infrared reflectance (A c ). 18 20 Note the difference 
in scale for the left and right hand axes: the anisotropy is 
actually decreasing with increased hole doping. The dashed 
lines are parabolic fits to the data, which are used to obtain 
interpolated values for doping levels not measured, and should 
be viewed as purely phenomenological. Lower panel: The 
ratio of the Lindemann number squared to the product of 
the in and out-of-plane penetration depths. This quantity, 
appearing on the right-hand side of Equation 2, controls the 
curvature of B m vs. T. 
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FIG. 5. The superconducting phase diagram of YBa2Cu30 !/ . 
H C 2(0) is suppressed in the same region that T c deviates from 
the parabolic form — a clear sign of the weakening of supercon- 
ductivity around 0.12 hole doping. The T c values are taken 
from Liang et al. 17 . 
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TABLE I. H C 2(0) and ex as obtained by fitting the vortex 
lattice melti ng curves to |Equation~2"l £o is calculated from 



H C 2(0) using Equation 1 The uncertainties come from the 
width of the resistive transition and the proximity of the low- 
est data point to T — K. The hole doping is obtained from 
the T c , using Figure 3 of Liang et al. 17 . 
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IV. DISCUSSION 

The suppression of T c in the underdoped region of the 
phase diagram was mapped in detail by Liang et al. 17 . 
In the same work, Liang et al. 17 correlated the c-axis 
lattice parameter with the hole doping of the copper- 
oxygen planes, showing a smooth evolution of hole doping 
with increased oxygen content. This demonstrates that 
the suppression of T c is not due to some peculiarity of the 
copper-oxygen chain doping mechanism in YBa2Cu30 y , 
but is in fact inherent to the electronic properties of the 
material. It was supposed that the suppression of T c 
may be due to a competition of superconductivity with 
stripe formation, as has been demonstrated explicitly in 
the lanthanum cuprates. 22 

The phase diagram in Fig ure 5| shows a clear correla- 
tion between the suppression of T c near 0.12 hole dop- 
ing and a suppression in the T = melting field and 
hence a suppression of H C 2(0). This further strengthens 
the case that the anomaly in T c is related to a weaken- 
ing of superconductivity. The corresponding maximum 
in coherence length — recall that £o ^ [H C 2{fy] — 
has also been seen in /iSR 23 and in the fluctuation- 
magnetoconductance 24 . 25 

Recent NMR 26 and x-ray diffraction 27 29 experiments 
have indicated the possibility of charge order in under- 



doped YBa2Cu30 y . In all three x-ray diffraction exper- 
iments, the charge order was seen to drop in intensity 
below T c . Additionally, Chang et al. 28 found that the in- 
tensity of the charge-order peaks could be increased with 
an applied magnetic field below T c . These experiments 
give further evidence for a close competition between su- 
perconductivity and the charge ordered state. This is in 
agreement with the minimum in H C 2(0) we observe near 
0.12 hole doping. 

V. CONCLUSION 

The onset of finite resistivity in a magnetic field coin- 
cides with the vortex melting transition in YBa 2 Cu3 0, / . 
This melting transition can be substantially below mean- 
field H C 2 at temperatures between K and T c . 12 Using 
a Lindemann criterion for melting produces good agree- 
ment between theory and experiment, with a Lindemann 
number cl between 0.3 and 0.4. These values are con- 
sistent with theoretical predictions, which vary between 
0.2 and 0.4 for highly anisotropic materials. 13 

Within the framework we used for flux-line-lattice 
melting, 7 12 13 B m is required to approach fJ,oH c2 as 
T — > 0. The agreement between our data and this theory 
suggests that )ioH C 2(0) = B m (0), in contrast to previous 
suggestions. 9 30 32 This means that the quantum oscilla- 
tions seen in underdoped YBa2Cu30 !/ would occur in a 
state free of vortices (superconducting fluctuations may 
still be present, 33 of course, as detected in the Nernst 
signal, 34 for example.) This absence of vortices is con- 
sistent with the lack of a field-dependent scattering term 
needed to fully describe the quantum oscillations. 35 

Below optimal doping, H C 2(0) is rapidly suppressed 
with decreasing hole doping, reaching a minimum of 24.5 
tesla at p = 0.116 holes. At lower hole doping H C 2(0) 
recovers — even as T c continues to decrease — indicating 
the presence of a phase that competes with supercon- 
ductivity, a phase which is strongest between 0.11 and 
0.13 holes. 
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